Cosmological perturbations of a perfect fluid and non-commuting variables 
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We describe the linear cosmological perturbations of a perfect fluid at the level of an action, 
providing thus an alternative to the standard approach based only on the equations of motion. This 
action is suited not only to perfect fluids with a barotropic equation of state, but also to those for 
which the pressure depends on two thermodynamical variables. By quantizing the system we find 
that 1) some perturbation fields exhibit a non-commutativity quite analogous to the one observed for 
a charged particle moving in a strong magnetic field, 2) local curvature and pressure perturbations 
cannot be measured simultaneously, 3) ghosts appear if the null energy condition is violated. 

PACS numbers: 04.20.-q, 98.80.Qc 



I. INTRODUCTION 

The theory of cosmological perturbations (TCP) for a 
perfect fluid has always been an important issue in cos- 
mology. It enables us to understand how small fluctua- 
tions seeded in the early universe eventually evolved into 
the present large scale structure. Also, TCP has been ex- 
tremely useful to put constraints on various cosmological 
models. 

TCP for a perfect fluid has been developed and stud- 
ied at the level of the basic equations of motion, i.e., 
the Einstein equations of general relativity (GR) and the 
energy-momentum conservation law [l], Q. Yet, TCP for 
a perfect fluid can be also studied at the level of the 
action. Although these two approaches are classically 
equivalent, the latter gives the following advantage. In 
TCP, one first has to perturb all the fields appearing in 
the equations of motion or in the action, such as the met- 
ric components and the energy density. However, as is 
well known, not all the perturbation fields are dynamical. 
Actually, GR with a perfect fluid has only one dynami- 
cal field for the scalar-type perturbation. But an identi- 
fication of this field as well as a derivation of its closed 
evolution equation by means of the equations of motion 
alone are not straightforward. 

The situation becomes worse when going to extended 
gravity models. For illustration, f(R, G) theories (R be- 
ing the Ricci scalar, and G the Gauss-Bonnet term) with 
a perfect fluid involve two dynamical fields for the scalar- 
tensor perturbation. In such theories, the usual approach 
based on the equations of motion requires a rather strong 
intuition because the closed evolution equations for those 
dynamical fields have to be extracted from rather compli- 
cated coupled differential equations. On the other hand, 
the action approach advocated in this paper allows a 
straightforward identification of the auxiliary fields just 
by checking the absence of any kinetic terms in the sec- 
ond order action. Once the auxiliary fields are found, 
they can easily be eliminated through their trivial equa- 
tions of motion. What is then left is an action contain- 
ing only the dynamical fields, from which we can derive 
the closed evolution equations. In [!, H[, we have explic- 



itly checked that the action approach indeed works for 
f{R, G) gravity models with no matter and with a scalar 
field, respectively. 

In this paper, we want to describe first-order TCP for 
GR with a perfect fluid at the level of the action, in a 
way consistent with the principles of thermodynamics. 
To this end, we use the action for a perfect fluid pro- 
posed by Schutz Q and do the quantization of the per- 
turbations, which might also be of some interest beyond a 
pure academic point of view. Indeed, the quantization of 
the background universe with a perfect fluid has been dis- 
cussed by many authors [fj|-[T3j. But here, we prove that 
quantizing the perturbation fields leads to non-standard 
commutation relations and, consequently, to unexpected 
effects upon the physical properties of any perfect fluid 
in quantum cosmology. 

A TCP action approach for fluids was first introduced 
in [l7j in the context of k-inflation. This approach was 
taken also in [l8[ to study non-linear cosmological per- 
turbations in the matter dominated universe. However, 
the fluid discussed there is the so-called scalar fluid whose 
energy-momentum tensor is completely written in terms 
of a scalar field and its derivative. By construction, the 
scalar fluid cannot have vector-type perturbations. Al- 
though there is an exact correspondence between a per- 
fect fluid and the scalar fluid for the scalar-type per- 
turbation at the linear order, it is no longer true for 
higher order perturbations because of the mixture of 
scalar and vector-type perturbations. On the other hand, 
the Schutz's action we will use here is for a perfect fluid. 
Therefore the action exactly describes the dynamics of 
a perfect fluid at any order. As far as we know, this is 
the first time TCP is fully developed within the Schutz's 
action. We believe our approach is suited for studying 
cosmology in extended gravity models. 

Before introducing the action for a perfect fluid, let us 
briefly review the thermodynamics needed to describe it. 
In this paper, we consider a "single" fluid, that is, a fluid 
whose thermodynamical quantities are completely deter- 
mined by only two variables, e.g. the chemical potential 
\i and the entropy per particle s [l4j]. In this sense one 
first needs to give two equations of state, n = n(fj,, s) and 
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T = T(p, s), where n is the number density and T is the 
temperature of the fluid. Using then the first law of ther- 
modynamics, dp = ndp — nTds, one obtains the pressure 
as p = p(p,s). Finally, the energy density is given by 
p = pn — p. This is enough to describe the system ther- 
modynamically. Single fluids also satisfy particle number 
conservation, namely N = nV is a constant. The second 
law of thermodynamics imposes d(Ns) — Nds > such 
that ds = at equilibrium. 

A single perfect fluid is also defined through its 
stress-energy tensor T^ v = (p + p)u^u v + pg^. In 
a Friedmann-Lemaitre-Robertson- Walker (FLRW) back- 
ground the conservation of energy-momentum, T^". v = 
0, implies that p + 3H(p + p) = or, equivalently, 
d(pV) + pdV = since V oc a 3 with a, the cosmolog- 
ical scale factor, and H = a/ a, the Hubble parameter. 
This, in turn, implies that diV = and ds = 0. In any 
FLRW universe we thus have na 3 = N = constant and 
s = 0. 



III. PERTURBATIONS 

Once and for all, we work within a spatially flat FLRW 
universe. At first order in perturbation theory we have 
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Note that Wi is gauge invariant since, following ref. [1| 
the perturbation fields transform respectively as 
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II. ACTION 

The action considered here has been introduced by 
Schutz [a and is defined as follows 



S = 



R 
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■p(fx,s) 



(1) 



Alternative functionals have been proposed, all being 
physically equivalent as shown in [15[. We chose the 
version (TTJ) as it was the most convenient for our pur- 
pose. The four-velocity of the perfect fluid is defined via 
potentials @: 



u v = — (d u £ + 6d u s + Ad v B) , 



(2) 



where £, 8, A and B are all scalar fields. The normaliza- 
tion for the four- velocity, u v u v = — 1, gives p in terms of 
the other fields. The fundamental fields over which the 
action (TTJ will be varied are g^ u , £, 8, s, A, and B. 
Having chosen the Lagrangian for gravity to be the 



one of GR. we recover G b 



inGT^ by varying with 



respect to the metric field. Besides the conservation of 
particle number and entropy already discussed^the other 
equations of motion derived from Eq. (TTJ) are [5(: 



u a d a 8 = T, u a d a A = 0, u a d a B = 0. 



(3) 



In a FLRW universe, ui = and uq 
solutions to Eq. ([3]) arc simply 



T such that the 



A = A(x) , B = B{x) , 



T{t')dt' + 8{x) . (4) 



There is a complete freedom for the functions A, B, and 8 
[2l| . any choice leading to the same physical background. 
We will take advantage of this freedom to simplify our 
study of the scalar and vector perturbations. 



under the gauge transformation x a — > x a + In 
Eq. (|6]) we have decomposed Wi into scalar (w s ) and 
divergence- less vector modes (v,i). So, in general Wi will 
generate both scalar and vector perturbations. However, 
we can efficiently use the freedom of choosing the time- 
independent background quantities A, B and 8 given in 
Eq. (0| to disentangle them. Any such choice does not fix 
a gauge as no conditions arc imposed on the perturbation 
fields themselves. 



A. Scalar type perturbations 



Let us simply consider the choice 



A = B = 8 = 0. 



(10) 



to remove the vector perturbations arising from Wi . Re- 
garding the metric, Sgoo and Sgoi are auxiliary fields such 
that the only scalar component which will be dynam- 
ical is the curvature perturbation (f> defined by Sgij = 
2a?4>8ij, with — > (f> — H£° under a gauge transforma- 
tion. 

We introduce the new quantity v = Si + 8(t)Ss such 
that Sui = di (v/p). Therefore, v represents the velocity 
perturbation of a perfect fluid. We then define two gauge 
invariant fields, <f> = <f> + Hv/p and 58 = 58 + Tv/p, 
to expand the action {T]) at second order, in a gauge- 
independent way: 
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The perturbation fields $ and 58 are related to the cur- 
vature and temperature, respectively. In the comoving 
gauge v = where a perfect fluid remains static, $ = (f> 
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and 59 = 59. The coefficients for the kinetic terms are 
given by [22j 



= P+P 2 _ P 

<S 9 7-7-9 ' ■ 



(12) 



whereas the remaining coefficients are 



C = 



H 



E = ^r, (13) 



and 



D = na 3 



(14) 



The general solution for 5s, 5A, and SB is their initial 
values since Eq. (|11[) forces them to be time-independent. 
As a consequence, the non trivial equations of motion are 
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a 3 Qs dt 
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These equations exactly coincide with those derived by 
perturbing the Einstein equations and the conservation 
law for the entropy, as it should be. In general, $ is 
sourced by 8s. For example, if the perfect fluid is an 
ideal non-relativistic gas characterized by T = mo), 
mo being the mass of the particles, then C ^ and we 
have to solve two coupled equations to know the time 
evolution of $ and 89. 

However, if T = f(s)fi, which is equivalent to hav- 
ing a barotropic equation of state p = p(p) [23| . then 
C = 0. (Note that both radiation and dust fulfill this 
condition, while a cosmological constant has vanishing 
Qs so that no contribution for perturbations arises, as is 
well known). In this case, the sign of Qs cannot be known 
from the usual approach based on the equations of motion 
alone. On the other hand, the action approach advertised 
here leads to an exact expression for Qs, which will be 
used to avoid ghost degrees of freedom when quantizing 
the perturbations. 

We also conclude that the fields $ completely decou- 
ples from 8s and propagates with a sound speed c s if 
C = and c 2 > 0. 



B. Vector type perturbations 

To arrive at the desired action via the shortest path, let 
us first assume that all the perturbation variables prop- 
agate only in one direction, say the z-direction. This 
should be allowed, as we know that perturbations with 
different wavenumber vectors do not mix in a FLRW uni- 
verse. Once we obtain the action for this particular mode, 
we can then easily infer the general action. 



The vector contributions come only from the compo- 
nent Ut of Wi defined in Eq. It is not easy to ex- 
tract Ui from this equation since the functions A, B and 
9 depend in general on the spatial coordinates. Yet, tak- 
ing again advantage of the freedom to select these back- 
ground functions, we can make the simplest choice that 
contains all the information needed for the vector modes, 
namely 



A = 9 = , B a = bi 



(17) 



where b = (b, 0, 0) is a constant vector orthogonal to the 
z-direction. With this assumption, we have w s = and 
Ui = bi SA(t, z) for Wi. 

Regarding the vector perturbation of the metric, we 
follow again ref. [l(| and denote Sgoi = aG;, and Sgij = 
a 2 (Cij+Cj t i), with transverse conditions G^i = Ci.i = 0, 
or, in our setup, G z = C' z = 0. Then, we impose the 
gauge condition 8B = 0. However, this condition alone 
does not completely fix the gauge, as only the component 
of £ l parallel to b gets frozen by Eq. Therefore we 
can still choose £ v such that C y = 0, and C = (7" is 
parallel to b. Finally, we find that the action for the 
vector perturbations is given by 

S v = J d 4 x[^^[(d z V x ) 2 + (d z V y f] +na 3 bSAC x 
+ na 2 bV x 8A + 2nGb 2 n 2 aSA 2 /#} , (18) 

where Vi = Gi — aGi is a gauge invariant field. This 
action can be immediately extended to the general case 
where the perturbation variables depend now on (x, y, z). 
In the gauge SB = 0, the result is given by 



d 4 o 
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(d J V l ){d ] V i ) + a 3 {p + p)C t Sm 



+ a 2 (p + p)Vi8ui - i a (p + p)8ui8ui] , (19) 



where we substituted Sut for biSA/fi. Variations with 
respect to Vi and Cj yield the following equations, 



AVi = 16irGa(p + p)5ui, 



dt 



[(p + p)a 3 8u % ] = 0, 



(20) 
(21) 



respectively. Again, these equations exactly coincide 
with those derived by perturbing the Einstein equations 
and the energy- momentum conservation law |16j |. This 
provides thus a cross-check that the calculations pre- 
sented here are correct. In fact, the main novelty in our 
approach is to be found when we quantize the system. 

To summarize section III, the known results on first- 
order TCP for a perfect fluid can be directly derived 
from variations of the classical action given in Eq. ((TJ). 
Note that a similar action approach has been already per- 
formed in [T3, [3- Yet, the system studied there cannot 
represent a perfect fluid. Indeed, as already mentioned in 
the introduction, the action proposed in [I?], EH is made 
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of a real scalar field. So, this system, by construction, 
cannot have vector perturbations, as the only new per- 
turbed field is the scalar one. Therefore the system stud- 
ied there is not a perfect fluid, otherwise a perfect fluid 
would have no vector perturbation. It can be thought of 
as a scalar fluid, but, once more, not as a perfect fluid. 
It is simply a different physical system whose squared 
sound speed c 2 is not equal to ft J p. 



IV. QUANTIZATION 

The most important advantage of the action approach 
proposed in this paper is, of course, that it allows us to 
quantize the system. Although the inhomogeneities of 
the present universe, such as the galaxy distribution, are 
clearly described by the classical theory, the quantiza- 
tion of a perfect fluid may have something to do with 
the early universe if the seeds for structure formation 
are provided by quantum fluctuations of fields generated 
during inflation. Yet, besides its practical utility, our ac- 
tion approach also opens new theoretical prospects, as 
discussed below. In the following, we will again treat the 
quantization for the scalar and vector type perturbations 
separately. 



A. Scalar type perturbations 

To quantize the scalar perturbations, let us first in- 
troduce the canonical field ip = \J a 3 Qs$. To avoid the 
appearance of a ghost, we assume that Qs is positive. Ac- 
cording to Eq. (fT2")l . this means that (p+p)/c 2 > 0. Such 
a constraint, together with the stability of the perturba- 
tions, c 2 > 0, lead to the null energy condition p + p > 0. 
Using the new variable ip, the action (jfip is rewritten as 



Ss 



d A x 



N - ■ D 

—(sess-Ssse) — Ss 



(22) 



where we have neglected SA and SB as they do not con- 
tribute to the Hamiltonian. The field ip has a canonical 
kinetic term, whereas the quadratic terms for 8 s and 80 
are at most linear in their time derivatives. Yet, it is 
known [l9| that a consistent quantization of such a sin- 
gular Lagrangian can be done provided one introduces 
the following equal-time commutation conditions, 



(t,x),7t(t,y)] = i8(x- y) , 



[Ss(t,x),S0(t,y)] 



N 



5(x-y) 



(23) 
(24) 



nian is given by 

H = I d 3 x 



2a- 



(25) 



One can easily check that the Heisenberg equations, with 
the help of the commutation relations, yield the same 
equations of motion as the classical ones derived from 
the variation of Eq. (|2"2"|) . 

The relation (|2"4")l shows that 5s and 80 become non- 
commuting variables at the quantum level. In Quan- 
tum Field Theory, different fields (i.e., different parti- 
cles) can be simultaneously observed at the same posi- 
tion. Here the perturbation fields related to the entropy 
and the temperature, to which we may individually at- 
tribute arbitrary numbers at the classical level, cannot be 
measured at the same space-time point. That this non- 
commutativity arises from the action of a perfect fluid is 
thus intriguing. 

We should concede that consequences directly linked to 
present observations are missing. However, at this level 
it is quite interesting to compare the action (|22p with 
the one of the Landau problem , an archetype of non- 
commutative geometry. Regarding 8s and 89, the action 
(|2"2"|) is essentially the same as the one for a charged parti- 
cle moving on a two-dimensional surface with a constant 
magnetic field background in the transverse direction: 



S= dt 



'•2 , -2^ B 



V ) - 2^ ±y ' 



yx) - V(x,y) 



(26) 



Within this analogy, the perturbation fields (8s, 80) cor- 
respond to the (x, y) space coordinates for the particle, 
and the number of particles N = na 3 plays the role of the 
constant magnetic field B. Interestingly enough, while 
the quite heuristic non-commutative relation [x,y] = 
— i/B in the Landau problem [f9| holds only in the ab- 
sence of the kinetic term in Eq. (|26[) . which is valid in 
the large magnetic field limit, the non-commutative rela- 
tion (f24| of a perfect fluid is exact for any finite number 
of particles. So, perfect fluids provide a nice example of 
non-commutativity between different fields. 

The other non-commutation relation (j2U)) leads also 
to an interesting physical consequence. By using once 
more the Einstein equations and the energy-momentum 
conservation law, we find that the pressure perturba- 
tion in the comoving gauge (y = 0) is given by 8p = 

— (p + p)4>/ H. Then, the commutator between <j) and 8p 
becomes 

[4>(t,£)Jp(t,jfl] =~ic 2 s H8{x-y)/a\ (27) 



All the other commutators are zero and 7r is the canonical 
conjugate momentum of The corresponding Hamilto- 



Consequcntly, local curvature and pressure perturbations 
cannot be measured simultaneously. 
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B. Vector type perturbations 

Time derivatives of Vi and Sui do not appear in the 
action (|1~9"1) . Therefore, those are auxiliary fields which 
can be eliminated through their equations of motion. The 
action (|19l) becomes then a functional which depends only 
on Cj. To make this action canonical, we introduce a new 
variable Fi(k,t) = V 'a 3 Q v (k,t)cf (k,t), where Cf(k,t) 
is the Fourier transform of Cj (x, t) and Qv is given by 



a 2 k 2 (p + p) 



k 2 + 16nGa 2 (p + p)' 



(28) 



To avoid the appearance of ghosts, Qv must be positive. 
So, as for the scalar modes we require p + p > 0, i.e. 
the null energy condition to hold. In terms of i 7 ), the 
canonical action in Fourier space is given by 

S v = [ dtd 3 k Uf^F, - W k F*F t ) , (29) 



with 

ml = -^log(a 3 Q v )-^ t loga 3 Q v 

Now the quantization is done by imposing the following 
canonical condition for i 7 ) and its conjugate momentum 



[Fi{t,k),if^{t,k')\ = i6(k-k') U 



The corresponding Hamiltonian is given by 



(31) 



H= I d 3 k Unfa + \m\F\F, ) , (32) 



and the evolution of the operators is given by the Heiscn- 
berg equation with the help of the commutation rela- 
tion (|3ip . The quantum version of Eq. ([2"U]) implies 
[Vi(t, x), 5uj(t, y)~\ = 0. Therefore, the gauge invariant 
metric perturbation and the vorticity of the perfect fluid 
can be measured at the same time, at the same position. 



As for the tensor perturbations, they come only from 
the metric perturbation. The action for the tensor per- 
turbations and its quantum aspects have been widely 
studied in the literature (e.g. [16(), mainly in connec- 
tion with the quantum generation during inflation. So, 
we do not discuss it any longer. 



V. CONCLUSIONS 

We have studied the theory of cosmological perturba- 
tions for a perfect fluid in GR at the action level. Starting 
from the action proposed by Schutz, we first reproduced 
the known results derived from the equations of motion 
alone. This enabled us to illustrate the advantage of our 
action approach at the classical level. Quantizing then 
the perturbation fields, we found that some of them do 
not commute, leading thus to a non-commutative field- 
geometry. In particular, we pointed out that a simulta- 
neous measurement of local curvature perturbations and 
pressure inhomogeneities is not allowed at the quantum 
level. Finally, we proved that both the null energy con- 
dition and a positive squared sound speed have to hold 
at all times in order to avoid ghost degrees of freedom. 

Another advantage of our action approach is that one 
can easily obtain the second order action depending only 
on the dynamical fields. Such an approach is thus suited 
to study cosmology in extended gravity models with more 
than one dynamical field. In particular, we expect that 
the approach presented here will be quite useful for the 
perturbation analysis of f(R,G) gravity models [20| . or 
for the treatment of non-gaussianities for the ent ropy and 
vector perturbations on perfect fluids following [1 81 ] . 



Acknowledgments 

We thank Sean Murray for helpful discussions. This 
work is supported by the Belgian Federal Office for Sci- 
entific, Technical and Cultural Affairs through the In- 
teruniversity Attraction Pole P6/11. 



[1] H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 78 
(1984) 1. 

[2] V. F. Mukhanov, H. A. Feldman and R. H. Branden- 

berger, Phys. Rept. 215, 203 (1992). 
[3] A. De Felice and T. Suyama, JCAP 0906, 034 (2009) 

[arXiv:0904.2092l [astro-ph.CO]]. 
[4] A. De Felice and T. Suy ama, Phys. Rev. D 80, 083523 

(2009) [arXiv:0907.5378l [astro-ph.CO]] . 
[5] B. F. Schutz, Phys. Rev. D 2, 2762 (1970) 
[6] P. Pedram and S. Jalalzadeh, Phys. Lett. B 659, 6 (2008) 

|arXiv:07U. 19961 [gr-qc]] . 
[7] P. Pedram, S. Jalalzadeh and S. S. Gousheh, Phys. Lett. 

B 655, 91 (2007) [arXiv:0708.4"l43l [gr-qc]]. 



[8] F. G. Alvarenga, J. C. Fabris, N. A. Lemos and 

G. A. Monerat, G en. Rel. Grav. 34, 651 (2002) 

[arXiv:gr-qc/01060"5l"] . 
[9] G. A. Monerat, E. V. Correa Silva, G. Oliveira-Neto, 

L. G. Ferreira Filho and N. A. Lemo s, Phys. Rev. D 73, 

044022 (2006) |arXiv:gr-qc/0508086| . 
[10] N. A. Lemos, J. Math. Phys. 37, 1449 (1996) 

[arXiv:gr-qc/9511082] . 
[11] V. D. Ivas hchuk and V. N. Melnik ov, Grav. Cosmol. 1, 

133 (1995) [arXiv:hep-th/9503223] . 
[12] P. Peter, E.' J. C. Pi nho and N. Pinto-Net o, Phys. Rev. 

D 73, 104017 (2006) |arXiv:gr-qc/0605060] 
[13] J. D. Brown, Phys. Rev. D 41, 1125 (1990). 



6 



[14] C. W. Misner, K. S. Thorne and J. A. Wheeler, San 

Francisco 1973, 1279p 
[15] B. F. Schutz and R. Sorkin, Annals Phys. 107 (1977) 1. 
[16] S. Weinberg, Oxford, UK: Oxford Univ. Pr. (2008) 593 

P 

[17] J. Garriga and V. F. Mukhanov, Phys. Lett. B 458, 219 
(1999) 

[18] L. Boubekeur, P. Creminelli, J. Norena and F. Vernizzi, 
JCAP 0808, 028 (2008) |arXiv:0806.10T6l [astro-ph]] . 

[19] R. Peierls, Z. Phys. 80, 763 (1933); L. Faddeev and 
R. Jackiw, Phys. Rev. Lett. 60, 1692 (1988); G. V. Dunne 



and R. Jac kiw, Nucl. Phys. Proc. S uppl. 33C, 114 (1993); 
R. Jackiw, |arXiv:hep-th/9306075| 
[20] A. De Felice, J. M. Gerard and T. Suyama, in prepara- 
tion. 

[21] Since u„ = (-1, 0), we also have that I = - f* fi(t')dt'+l, 
and Vt = -AVB, which implies that VA x VB = 0. 

[22] For c\ we used the fact that p = (dp/dp) s p+ (dp/ds) s. 

[23] In this case, we obtain (d/x/ds) = T such that 
(dp/ds) p = n[(d f i/ds) p - T] = 0. 



